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Without invoking a perturbative expansion, we define the cosmological curvature perturbation, 
and consider its behaviour assuming that the universe is smooth over a sufficiently large comoving 
scale. The equations are simple, resembling closely the first-order equations, and they lead to 
results which generalise those already proven in linear perturbation theory and (in part) in second- 
order perturbation theory. In particular, the curvature perturbation is conserved provided that the 
pressure is a unique function of the energy density. 
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I. INTRODUCTION 



In this note we define and study a non-linear generalisation of the linear cosmological 'curvature perturbation' 
(f-) ' which is commonly denoted by (. Just as in linear theory, we find that the curvature perturbation is conserved during 
any era when the locally-defined pressure is a unique function of the locally-defined energy density, but is otherwise 
C*~) time-dependent. 

Before describing the calculation we mention the motivation which is actually two-fold. The first motivation 
concerns the comparison of theory with observation. Observation has now established the existence of a nearly 
| Gaussian and scale-invariant spatial curvature perturbation, present already before comoving cosmological scales 
come inside the Hubble distance (enter the horizon). The curvature perturbation presumably is generated somehow 
^-H ' by the perturbation of a scalar field, that perturbation in turn being created on each scale at horizon exit during 
inflation. Several mechanisms for achieving this have been proposed. In all of them it is supposed that the curvature 
perturbation, defined on suitable spacetime slices, is conserved (starting from the time when the creation mechanism 
ceases) until the approach of horizon entry. The conservation is supposed to hold by virtue of the condition, either 
Q H - proven or just assumed, that the pressure is a unique function of the energy density. This is the 'adiabatic pressure' 
■ ' condition, that in first order perturbation theory leads to 5P/5p = P/f>- 

The conservation of the curvature perturbation for adiabatic pressure has been demonstrated at linear order for an 
appropriate slicing of spacetime. That slicing can be taken to be the either the 'comoving' one (the one orthogonal to 
comoving worldlines) or the uniform-density one or the uniform-proper-expansion (uniform-Hubble) one; it does not 
matter which because in linear theory the three are known to coincide in the large-scale limit P, Q . 

These results are all that one needs at the moment, but with better data in the future it will be necessary to go to 
second order. This is because primordial non-Gaussianity, which may well be big enough to be detected, is typically 
generated only at the second-order level. The study of second-order gauge invariant perturbations started only quite 
recently 0,0] and has been hampered by it's sheer complexity. The prospect of observing non-Gaussianity in the near 
future has led to renewed activity in this area [U, 0, 01 H, 13 but the second-order approach has not yet been carried 
through to the point where it sufficiently generalises the first-order large-scale results. 

As serious as its incompleteness is the fact that even on very large scales the second order calculation is complicated, 
losing as a result the physical transparency of the first-order calculation. This brings us to our first motivation; we 
will provide a non-linear understanding of the large-scale situation (applying therefore to all orders in the perturbative 
expansion) which follows the same lines as the first-order discussion and is not much more complicated. We will show 
explicitly how to recover powerful second-order results from this approach, without going through the complicated 
second-order formalism. 

Our second motivation is more philosophical, concerning the possibility that inflation lasted for an exponentially 
large number of e-folds. In that case the expanding Universe which we observe is part of a huge region, which at 
the classical level was homogeneous before it left the horizon. When comparing the theory with observation, very 
distant parts of this huge region are irrelevant, because one can and should formulate the theory within a comoving 
box whose present size is just a few powers of ten bigger than the Hubble distance. Linear theory with the small 
second-order correction is then valid within the box, which is all that one needs for practical purposes. Still, it is of 
interest to understand the nature of the universe in very distant regions. In this connection it is important to recognise 
that a Gaussian 'perturbation' with a scale-invariant spectrum actually becomes indefinitely large in an indefinitely 
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large region. This means that one should try to understand the cosmological perturbations without invoking the 
perturbative expansion, which is precisely what we do here. 

We work with the ADM formalism jljj]- Historically, the ADM-formalism was employed to study non-linearities 
already by Bardeen in albeit in the context of perturbation theory. Non-perturbative studies of fluctuations 
followed in ( see also |l2l In the relation between PPN-formalism and linear perturbation theory was 

established using the ADM-formalism and in 0] the formation of primordial black holes was studied numerically. 

In addition, we should mention the existence of related work done by the Russian school, which is based on 
construction of solutions of Einstein equations having locally Friedmann-like behaviour near the singularity, the so- 
called " quasi-isotropic" approach. It was pioneered by Lifshits and Khalatnikov in 1963 [l6|. generalised to the de 
Sitter case and has been elaborated recently 3 El • 

This paper is organised as follows. In Section^ we define the metric of our spacetime, specify the assumptions 
for its validity, and describe the basic geometrical properties under our assumptions. In Section II I II we show the 
conservation of a non-linear generalisation of the curvature perturbation on very large scales. In Section IIV | we 
express our conserved quantity in a slice-independent manner and relate it to the AN formula derived in [2i| (see 
also We conclude in the final section. In Appendix A we give the components of the relevant tensors and in 

Appendix B we give the spatial traceless components of the Einstein equations. 

II. GEOMETRY AND ENERGY CONSERVATION 
A. Metric 

We use the standard (3 + l)-decomposition of the metric (ADM formalism), which applies to any smooth space- 
time [13: 

ds 2 = -M 2 dt 2 + jij {dx l + f3 l dt)(dx J + ftdt) , (1) 

where Af is the lapse function, j3 l the shift vector, and 7y the spatial three metric. (Greek indices will take the values 
fj,, v = 0,1,2,3, Latin indices i,j = 1,2,3. The spatial indices are to be raised or lowered by r ) lJ or 7y.) In this 
(3 + l)-decomposition, the unit time-like vector normal to the x° = t =const. hypersurface n M has the components, 



1 (3 l 

a7'~a7 



(2) 



n M = [-Af, 0] , n" 
We write the 3-metric, 7^, as a product of two terms, 

lij = e 2a 7y , (3) 

where a and 7^ are functions of the spacetime coordinates (t,x l ), and det[7y] = 1. Because of the latter condition, 
the first factor is a locally-defined scale factor which we denote by a, 

e a = ~a. (4) 

We are interested though in the inhomogeneity of a, and so we factor out from e Q a global scale factor a(t) and a 
perturbation ip, 

e a = a(t)e^'> xi ) . (5) 

We assume that ip vanishes somewhere in the observable Universe (say at our location). This makes a(t) the scale 
factor for our part of the universe, and ensures that ip is a small perturbation throughout the observable Universe. 
(If we are surrounded by a super-large region within which ip becomes large, then the appropriate scale factor for a 
generic observer should be a evaluated at their own location, making again ip small in their vicinity.) 
In a similar way the matrix 7^ can be factored, 

7 = Ie H , (6) 

where / is the unit matrix. The condition detifj) = 1 ensures that the matrix H is traceless, which follows from the 
relation det(exp(M)) = exp(tr(M)), valid for any symmetric matrix M. In our part of the universe, with coordinates 
corresponding to the usual gauges, Hij is a small perturbation. 
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B. Gradient expansion 

Cosmological perturbation theory expands the exact equations in powers of the perturbations, keeping only terms 
of a finite order. In particular, first-order perturbation theory linearises the exact equations. We will not use this 
perturbative approach, but instead will use the gradient expansion method 0,0, H^l, which is an expansion in 
the spatial gradient of these inhomogeneities. To be precise, we focus on some fixed time, and multiply each spatial 
gradient di by a fictitious parameter e, and expand the exact equations as a power series in e. Then we keep only the 
zero- and first-order terms and finally set e = 1. 

In the perturbative approach the fictitious parameter e would be multiplying the perturbations. Working to linear 
order in the e of the gradient expansion obviously reproduces that subset of the linear perturbation theory expressions 
which can be derived by considering only those equations which have at least one spatial gradient acting on each 
perturbation. This is why our results will closely resemble those of linear perturbation theory. 

The gradient expansion is useful when every quantity can be assumed to be smooth on some sufficiently large 
scale with coordinate size If we are to model our actual universe by this smoothed universe, it is necessary to 
implement a smoothing procedure at the level of the field equations, either of Einstein gravity or of any alternative 
theory. However, smoothing the smaller scale inhomogeneities is a delicate issue, which is beyond the scope of the 
present paper. Here, we simply assume that there exists some kind of smoothing that can give a good approximation 
to the actual universe on coordinate scales greater than fc -1 . 1 Focusing on the observable Universe, this corresponds 
to a comoving smoothing scale of physical size a(t)/k. Then, instead of introducing e as a formal multiplier of the 
spatial gradients, we can make the identification 

e=k/aH, (7) 

where H = a/a is the Hubble parameter. 2 At a fixed time the limit e — > corresponds to k — ► 0. In Hubble units, 
the typical value of the gradient of a quantity / will be ef. 

Our key physical assumption is that in the limit e — > 0, corresponding to a sufficiently large smoothing scale, 
the universe becomes locally homogeneous and isotropic (a FLRW universe). (By 'locally' we mean that a region 
significantly smaller than the smoothing scale, but larger than the Hubble scale, is to be considered.) The Hubble 
distance is the only geometric scale in the unperturbed universe. In the perturbed Universe there is in addition the 
scale 1/k under consideration, and possibly other scales provided by the stress-energy tensor. Unless one of the latter 
scales is bigger than 1/k, local homogeneity and isotropy will be a good approximation achieved throughout the entire 
super-horizon era k -C aH, and the results of this paper will be valid throughout that era. Because cosmological scales 
are so large, one expects this 'separate universe' hypothesis |24l l25l |2(| to be valid for them, ensuring the maximum 
regime of applicability for our results. 

An immediate consequence of our assumption is that the locally measurable parts of the metric should reduce to 
those of the FLRW. Thus there exists an appropriate set of coordinates with which the metric of any local region can 
be written 

ds 2 = -dt 2 + a 2 {t)5 lJ dx l dx 1 . (8) 

(We took this metric to be spatially flat, which is the expectation from inflation and agrees with observation; a small 
homogeneous curvature would make no difference.) 

Let us see what this implies for the metric components. In the limit e — > 0, the above local metric should be globally 
valid. This implies that the metric component Pi vanishes in this limit, (3i — 0(e). 3 It may be noted, however, that 
this is not really a necessary condition but rather a matter of choice of coordinates for convenience. 

What about the quantity 7^? A homogeneous time- independent 7^ can be locally transformed away by choice 
of the spatial coordinates, but a homogeneous time-dependent 7^ is forbidden because it would not correspond to a 
FLRW universe. We therefore require 7 = O(e). In Appendix 151 though, we show that 7 decays like a~ 3 in Einstein 
gravity if it is really linear in e. Taking the usual view that decaying perturbations are to be ignored, we conclude 
that 7 will be of second order in e. 



1 In linear perturbation theory it is appropriate to use a Fourier expansion and then smoothing corresponds to dropping Fourier components 
with wavenumber bigger than fc, but we have no use here for the Fourier expansion. Still, it is useful to keep the case of the Fourier 
expansion in mind. 

2 If we are surrounded by a super-large region, a should be replaced by the local scale factor a(t, x 1 ), and H by the local Hubble parameter 
H(t,x r ) that we define later. 

3 We adopt the traditional mathematics notation l27l . according to which / = 0(e n ) means that / falls like t n or faster. 



4 



The conditions on the metric components are therefore 

ft = 
% = 



0(e), 
0(e 2 ). 



(9) 
(10) 



There is no requirement on ip and Af since they are not locally observable. We note that in alternative theories 
of gravity, the assumption 7 = 0(e 2 ) may not be as natural as in the Einstein case. Nevertheless, we assume this 
condition. In other words, we implicitly focus on a class of gravitational theories in which the condition 7 = 0(e 2 ) is 
consistent with the field equations. 

In view of Eq. 0, the line element simplifies, giving 



ds 2 = -M 2 dt 2 + 2f3 l dx l dt + j.jdx'dx 3 
C. Energy conservation 



(11) 



By virtue of the separate universe assumption, the energy momentum tensor will have the perfect fluid form 

T MV = (P + P) u M u v + g^P , (12) 

where p — p(x^) is the energy density and P = P(x M ) is the pressure. 

First let us choose the spatial coordinates that comove with the fluid. That is, the threading of the spatial 
coordinates such that the threads x % =constant coincide with the integral curves of the 4-velocity (the comoving 
worldlines). Hence, 



dx 1 
~~dt 



The components of the 4-velocity in these coordinates are 

1 



.0 



= 0. 



+ 0(e 2 ), 



(13) 



ft 



^JM 2 ^WWk_ 

The expansion of u M in the comoving coordinates, v l = 0, is given by 



-A/", 



Af 



+ 0(e 2 ). 



= V M u" 



-.d u 



3^0 (AA 



,3a„,0\ 



Aft 



,3a 



Afe 3a \ ^M 2 - Wi 



(14) 



(15) 



Note that 7^ does not appear in the above expression because det jij = 1. The relation between the coordinate time 
x° = t and the proper time r along u M is 



dt 



= u° = 



1 



yJM 2 - Pfc 



The energy conservation equation, 



dr 



P + (P + P) ' 



= 0, 



reduces therefore to 



where 



dr 



P+(P + P) 



3d 

77 



= p + 3{p + P)a + 0(e 2 ) = 0, 



0(e 2 



(16) 



(17) 



(18) 



(19) 
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It is important to note that the expansion of the hypersurface normal n M is given by 

9 n ee V,n" = (e 3 «/3*) . (20) 

Thus 8 and 0„ are equal to each other at linear order in e. Note that the above argument uses only the energy 
conservation law, hence applies to any gravitational theory as long as the energy conservation law holds. 

Here let us point out a couple of immediate but important implications of the above. The first point is that the 
equivalence of 9 and 9 n for the comoving threading readily implies the equivalence of 9 and 9 n for any choice of 
threading for which f3 l = O(e). This is because the change of the threading affects the numerical value of 9 n at 
a given world point (t,x l ) only by terms of 0(e 2 ) as is clear from Eq. (|20|l . The second point is that, the above 
argument applies not only to the total fluid but also to any sub-component of the fluid, provided that it does not 
exchange energy with the rest, and that the comoving threading with respect to that component satisfies the condition 
j3 l = O(e), that is, if the 3-velocity v l remains 0(e) for any threading with j3 % = 0(e). 

Once we have the equivalence of 9 and 9 n , it is convenient to introduce the notion of a local 'Hubble parameter' 
3H = 9„, 

H = \9 n = £(£-Hfr) +0 (e>). (21) 
As is shown in Appendix^ we then recover a local Friedmann equation once we appeal to the Einstein equations. 



III. THE EVOLUTION OF THE CURVATURE PERTURBATION 



Now let us investigate the evolution of the curvature perturbation ip. So far, we have not specified the choice of the 
time-slicing. Below we shall consider some typical choices of the time-slicing separately. 



A. Uniform-density slicing 

In this subsection we consider the uniform-density slicing, denoting ip on this slicing by — £. We shall need only the 
condition (3 — 0(e), not the other condition 7 = 0(e 2 ). 

Following the linear treatment of [23 , we avoid in this subsection the assumption of Einstein gravity. Instead we 
just consider some energy-momentum tensor T^, which satisfies V ' V T^ V = corresponding to energy-momentum 
conservation. This is useful in two ways. First, it allows us to deal if desired with just one component of the cosmic 
fluid instead of the total. Second, our results will apply to the case (arising for instance in RSII cosmology [2£|) that 
Einstein gravity is actually modified. 

Throughout this paper we are working to first order in e. We take the anisotropic stress of the fluid to be negligible 
at that order. In other words, the anisotropic stress is supposed to be of second order in e. (This can be verified in 
specific cases, in particular if the fluid consists of a gas and/or scalar fields.) Then, because there is a unique local 
expansion rate H to linear order in e, the local energy conservation equation (|18fl has the unperturbed form to this 
order, 

±p=-3H(p + P) + 0(e 2 ). (22) 

Multiplying each side by Af this becomes 

-+i> = -\^ + 0(e 2 ). (23) 
a 3 p + P 

At each point this equation is valid independently of the slicing. Now let us go to the uniform-density time-slicing, 
and denote ip on this slicing by — £. If, to first order in e, P is a unique function of p (the 'adiabatic pressure' condition), 
then Eq. (|23(l shows that tjj is spatially homogeneous to first order. Since ip is supposed to vanish at say our position, 
though the position can be chosen arbitrarily, this means that ip on uniform-density slices is time-independent to first 
order, 

-^ = C = 0(e 2 ). (24) 

At this stage, as noted before, £ can refer to the total cosmic fluid, or to a single component which does not exchange 
energy with the remainder. 
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B. The comoving and uniform-Hubble slicings 



Now we invoke Einstein gravity, and show that the comoving and uniform-Hubble slicings coincide to first order in 
e with the uniform-density slicing. By 'comoving slicing' we mean the one orthogonal to the comoving worldlines (it 
should be perhaps called the velocity-orthogonal slicing from the general relativity point of view, but the terminology 
'comoving slicing' is usual in cosmology). As we have already decided to use the comoving worldlines as the threading 
this fixes the gauge completely, and we call it the comoving gauge. 

By invoking the comoving slicing we are setting the vorticity of the fluid flow equal to zero, since the slicing exists if 
and only if that is the case. This is reasonable because vorticity is not generated from the vacuum fluctuation during 
inflation. Also, for a perfect fluid, there is a vorticity conservation law in arbitrary spacetimc 29], which states that 
the magnitude of the vorticity vector is inversely proportional to S'exp^ dP/(p + P)] along each fluid line, where S 
is the cross-sectional area of a congruence of the fluid orthogonal to the vorticity vector. Thus, the vorticity would 
become negligible a few Hubble times after horizon exit on each scale, even if it were somehow generated during 
inflation. 

The Einstein equations are 

G>„ = SnGT^v , (25) 

where refers now to the total fluid, and G M „ and G are the Einstein tensor and Newton's constant, respectively. 
We consider its components in Appendix A3. The (0, z)-component gives in the comoving gauge 

diH = 0(e 3 ) , (26) 
and the (0, 0)-component gives the local Friedmann equation 

H2 = ^ p + (e 2 ), (27) 

leading to dip = 0(e 3 ). This makes the typical magnitudes of Sp and 8H go like e 2 . In other words, the uniform- 
density, uniform Hubble and comoving slices coincide to linear order in e. These results are the same as those of linear 
perturbation theory. 

Knowing that H and p are both spatially homogeneous, we can learn about the lapse function by writing the energy 
conservation equation Eq. (|22|l in terms of coordinate time, 

l p =-3H(p + P) + 0(e 2 ). (28) 

Since p and H are spatially uniform we learn that J\f is of the form 

^=*TTpb) +0< « 2) ' (29) 

where A(t) can be any function which makes J\f positive definite. If the pressure is adiabatic, J\f is independent of 
position and we can choose Af = 1. In that case, the comoving gauge is also a synchronous gauge to first order in e. 
If the pressure is not adiabatic we can write on the uniform energy density slicing 

n= (so) 

This is the non-linear generalisation of the known result in first-order perturbation theory, 

»-i-£p. (=") 

where SP is the pressure perturbation on uniform-density slices (the 'non-adiabatic' pressure perturbation). 4 

As shown in Appendix iBl the traceless part of the spatial components of the Einstein equations is 0(e 2 ), while the 
trace part gives no additional information. 



4 As we were preparing the present paper, a related one appeared |30H claiming that the spatial variation of the coordinate expansion rate 
in the non-adiabatic case could change the presently-accepted predictions for observable quantities. We disagree with this conclusion, 
since that variation is already implicitly present in any correct formulation of cosmological perturbation theory. 
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IV. GAUGE TRANSFORMATIONS AND AN FORMULA 



As we are working to first order in e, the threading is unique in the sense that all the threadings are equivalent to 
the comoving threading as discussed in Section III CI leaving only the slicing to be fixed. In this section we consider 
the effect of a change of slicing on the curvature perturbation. 

Let us define the number of e-foldings of expansion along an integral curve of the 4- velocity (a comoving worldline): 



N{t2,h]X*) = - 



eNdt = -- 

3 



dt 



P 



(32) 



where, for definiteness, we have chosen the spatial coordinates {x 1 } to be comoving with the fluid. The essential point 
to be kept in mind is that this definition is purely geometrical, independent of the gravitational theory one has in 
mind, and applies to any choice of time-slicing. 
From Eq. (|2*T1) we find 



i>{t 2 , x l ) - Mux*) = Nfatux*) - In 



o(t a ) 



(33) 



Thus we have the very general result that the change in tjj, going from one slice to another, is equal to the difference 
between the actual number of e-foldings and the background value No(t 2 ,ti) = ln[a(t 2 ) / a(ti)]. One immediate 
consequence of this is that the number of e-foldings between two time slices will be equal to the background value, if 
we choose the 'fiat slicing' on which ip — 0. (This slicing is of course truly flat only if -yy = &y.) Thus the flat slicing 
is one of the uniform integrated expansion slicings |26j . 

Consider now two different time-slicings, say slicings A and B, which coincide at t = t\ for a given spatial point 
x l of our interest (i.e., the 3-surfaces £^(ti) and Es(ti) are tangent to each other at x l ). Then the difference in the 
time-slicing at some other time t — t 2 can be described by the difference in the number of e-foldings. From Eq. I|33|l . 
we have 



«M*2, x*) - ip B (t 2 ,x l ) 



JVkfta.tija;*) - N B (t 2 ,tv. 
&N AB {t 2 ,x % ), 



(34) 



where the indices A and B denote the slices A and B, respectively, on which the quantities are to be evaluated. 
As discussed in Section |V| this generalises, for large scales only, the known result of the first- and second-order 
perturbation theory 0, • 

Now let us choose the slicing A to be such that it starts on a flat slice at t = t\ and ends on a uniform-density slice 
at t = t 2 , and take B to be the flat slicing all the time from t = t\ to t = t 2 . Then applying Eq. Q34JI to this case, we 
have 



ip A (t 2 ,x l ) = N A (t 2 ,t 1 ;x i )-N (t 2 ,t 1 ) = ANpit^x 1 ). 



(35) 



where ANp(t 2l t\; x l ) is the difference in the number of e-foldings (from t = t\ to t = t 2 ) between the uniform-density 
slicing and the flat slicing. This is a non-linear version of the AiV formula that generalises the first-order result of 
Sasaki and Stewart pOj . 

Now we specialise to the case P = P(p). In this case, Eq. I|33l) reduces to 



^(ta.a^-Vfa.s*) 



111 



~a(t 2 y 




r-p(t 2 .X l ) 


dp 




3 


/(o(ii,K 4 ) 


p + p 



Thus, there is a conserved quantity, which is independent of the choice of time-slicing, given by 



CO**) 



p(t,x') 

pit) p 



dp 



P 



(36) 



(37) 



In the limit of linear theory, this reduces to the conserved curvature perturbation in the uniform-density, uniform- 
Hubble, or the comoving slicing, 



((x i )=K c (x l ) = ^(t,x l ) + 



3(P + P) 



(38) 



Finally, we mention that the generalisation of all the above results to the case of arbitrary threading, not restricted 
by the condition f3 l = O(e), is formally trivial. Let us denote the general spatial coordinates by {X 1 }. They are related 
to the comoving coordinates {x 1 } by a set of coordinate transformations, x % = F l (t, X 1 ). Then all the equations above 
are valid for an arbitrary choice of threading by simply replacing the arguments x % of all the functions by F l (t,X l ). 
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V. APPLICATIONS TO COSMOLOGICAL PERTURBATION THEORY 

In this section we make contact with cosmological perturbation theory, showing how our results both reproduce 
and extend known second-order results. 

Let us begin with the first-order case. To first order in the perturbations the spatial metric becomes 5 

QH = a 2 (t) [Sij (1 + 2^j) + Hij] , (39) 

In this context it is known (see for instance [2(|) that (always referring to the super-horizon regime) three slices coin- 
cide: uniform density, uniform proper expansion (uniform Hubble, the expansion being independent of the threading), 
and comoving slice. Also it is known that the perturbation ip is independent of the threading. Finally, defining the 
curvature perturbation £ as the value of — ip on this 'triple coincidence' slicing, it is known that C is conserved as long 
as pressure is a unique function of energy density. 

Going to second order (with an appropriate definition of ip) onl y so me of these statements have been verified. In 
particular, the conservation of C, was shown by Sasaki and Shibata |l5j] in non-linear theory if P/p is constant and by 
Malik and Wands [3 in second order perturbation theory setting 7.^ = Sij (i.e. ignoring the tensor), and by Salopek 
and Bond during single-component inflation. The main point of our paper has been to show that all of them are in 
fact valid (with again an appropriate definition of ip). 

Our finding is important because it ensures that the non-gaussianity of the curvature perturbation can be calculated 
once and for all at the time of its creation, remainin g th ereafter constant until horizon entry [jj^] . This was the implicit 
assumption made by Maldacena @ (see also |83ll84l IsH Isfij): he calculated the non-gaussianity of the curvature 
perturbation (to be precise, its bispectrum) a few Hubble times after horizon exit in a single-field model, defining 
it on the comoving slicing. It was also that of |37|. who calculated the curvature perturbation just before curvaton 
decay, now on the uniform density slicing, and that of |38| who calculated it at the end of inflation in a two-component 
inflation model with a straight inflaton trajectory (again on the uniform-density slicing). 

The gauge transformations and gauge- invariant expressions that we derived in the last section reproduce the second- 
order results. To see this, let us first consider two definitions of the curvature perturbation in the literature. 6 In all 
cases we shall employ the notation that a generic perturbation g is split into first- and second-order parts according 
to 

9 = 9i + ^92 ■ (40) 

One definition, used by Maldacena (introduced by Salopek and Bond in 11]) to calculate the non-Gaussianity gener- 
ated by single-field inflation, coincides with our definition of tp, 



2a = a 2 (t)e 2C = a 2 {t) (l + 2( + 2( 2 ) . (41) 



The other generalisation, employed by Malik and Wands (based on is different and we denote it by (mw; 
it is 

e 2a = a 2 (t)(l + 2Cmw) , (42) 

so that 

Cmw = C + C 2 , (43) 

or equivalently 

C2mw = C2 + 2(G) 2 . (44) 
Evaluated to first order in the perturbations, the gauge transformation Eq. (|34|l reduces to the known result 

ipA ipB — HAt , (45) 



5 The notation tp is that of |3ll : Kodama and Sasaki Jjj denote the same quantity by TV). 

6 Another definition is discussed elsewhere \,\A . 
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where At is the time displacement between the slices. If one of the slicings has uniform density and the other is flat 
this gives a gauge-invariant definition of C, 

-C = ^ + ^-Sp, (46) 
Po 

where a prime denotes differentiation with respect to conformal time and TL = a' /a, and the right hand side is 
evaluated on a generic slicing. 

To illustrate how things work at second order, we will jus t show how Eq. (|34() reproduces the known second-order 
result 7]. 7 Following along the lines of Lyth and Wands |2g, we expand the integrated expansion to second order in 
a power series expansion centred on the flat slicing of Eq. (|34|) , 

Using Eqs. and g7} this gives the perturbed expansion at second order, 

SN 2 = %5p 2 2^5p[5 Pl + (n P j H') feY , (48) 
Po Po V Po / V Po J 

where the right hand side is evaluated on flat slices. The expression given in Eq. 148JI is then related to C2mw, if this 
is also evaluated on flat slices, by C2mw = —SN 2 + 2Ci, which coincides with Eq. (|4^|l using Eq. J57JI, i.e. (2 = —SN 2 . 



VI. CONCLUSIONS 



Our central result is the existence of a conserved non-perturbative quantity £, corresponding to the scalar curvature 
perturbation in the perturbative case, which may be defined on the uniform-density, uniform-expansion or comoving 
slices since these coincide in the large scale limit. Locally, this statement follows from the equations of the coordinate- 
free approach as given for instance in |39j and reviewed in |4(ij|. 8 We have here preferred to employ the usual coordinate 
approach. This is because the coordinate approach is used in practise in almost all treatments of the evolution of 
perturbations during and after horizon entry, owing to the relative ease with which it handles the effect of particle 
collisions and free-streaming. Using the coordinate approach has allowed us to make contact with existing second 
order perturbation calculations, clarifying some previously mysterious connections between those of different authors. 
Also, on very, very large scales, it makes contact with the idea that the 'perturbations' are actually supposed to 
be random fields which can become arbitrarily large in an arbitrarily large region, thus providing in some sense a 
generalisation of the stochastic description of scalar field evolution during inflation. Finally, the gradient expansion 
can be most easily implemented in the coordinate approach, in particular if we want to go to the next order in the 
gradient expansion. 
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7 That this happens is stated without proof in , and for the special case of adiabatic pressure in 12(1 . 

8 Wh ile finishing version 2 of the present paper, a related work on the evolution of perturbations appeared, using the covariant approach 

El. 
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The metric tensor is given by 



APPENDIX A: TENSOR COMPONENTS 



The metric tensor 



.goo = -A/" 2 + ft ft , g Ql = g l0 = ft 
9ij = lij = e Ty > 



and 



where detpyy] = 1 and ft — ft 



..on 



ft ft 



<? Oi = <7 i0 



— ^3 

,2a ' 



TV 2 ' 
ft ft 
In 



3-geometry 



(Al) 
(A2) 



(A3) 
(A4) 



The unit timelike vector normal to the hypersurface t =const. is given by 



= [-AT, 0,0,0] , 



1 ft 1 



This gives for the extrinsic curvature tensor [42 



1 

w 



—fifty + % + Put 



which can be expressed as 



K i:j = -Y^ij + Aj , 



where 7 y ^4y = and 6 n is the expansion of n p 



On = V M n" . 



On large scales we have Aij — 0(e) by virtue of our assumptions, and Eq. I|A8J) reduces to 



+vu: + o( e ^). 



(A5) 

(A6) 
(A7) 

(A8) 
(A9) 
(A10) 



In the case of a conformally flat 3-geometry, i.e., for 7^ = 5ij, the intrinsic curvature on spatial 3-hypersurfaces is 
expressed 



(All) 



Einstein tensor 



Here we give the components of the Einstein tensor to first order in e, valid in an arbitrary gauge provided that the 
metric satisfies the conditions © and l]10p. 

First let us consider the case of ft — 0. Then, the time derivative d/dt is along the normal vector n u = (—Af, 0). 
Then, for the (0, 0) and (0, i) components of the Einstein tensor, we can apply the standard (3 + l)-decomposition of 
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the Ricci curvature tensor (i.e. the Gauss-Codacci equations), which are essentially the Hamiltonian and momentum 
constraint equations. They are given, for example, by Eqs. (2.20) and (2.21) of In our notation, they are 

G% = -\ ( (3) i? + \k 2 - A ij Aij \ , (A12) 
G°j = D i A i j -^DjK , (A13) 



where Di is the covariant derivative with respect to the metric 7^ , and the extrinsic curvature Kjj is given by Eq. (|A8I) 
with 6 n = —K. 

Since Ay = 0(e 2 ), and ^>R involves at least second derivatives of the metric tensor, these expressions reduce to 

G°o = - l -K 2 + 0{e 2 ) 1 (A14) 

G° 3 = -^l Dj K + 0(e 3 ), (A15) 

where 

K = -^(^ + ^j+0(e 2 ). (A16) 

For the components, they can be decomposed into the trace and traceless part. Evaluation of the traceless 
part is slightly involved. We write down the corresponding components of the Einstein equations in Appendix [B] hi 
which it is shown under our assumptions that they are 0(e 2 ) and hence can be neglected. The trace part reduces to 

G\ = jfdtK - K 2 . (A17) 



Therefore 



G l 3 = \ (jjd t K - K 2 ^j 5) + 0{e 2 ) . (A18) 

It is easy to see that inclusion of /3j = O(e) does not change these results at all. The above results can be regarded 
as the (3 + l)-decomposition with the hypersurface normal vector n" = (1/Af, 0). Thus, for a general choice of the 
spatial coordinates, we just have to perform the following replacements: 

G% -> -»„(?",»" = G% + G°# , (A19) 
AfG°j -n lt G' i v h v j = AfG°j , (A20) 
G*j -» h* ^G" vh 1 - 'j = G° j (3 i + G i j . (A21) 

So, the difference is of the form, G°j(3 l , which is 0(e 2 ). 

To summarise, the components of the Einstein tensor on large scales are given by 

G% = - l -K 2 + 0( ( 2 ) 1 (A22) 

G°j = -^\D 3 K + 0{e 2 ), (A23) 

G l 3 = I (jjdtK - K 2 ^j 8) + 0(e 2 ) , (A24) 

where K is given by Eq. JM6j|. 

Although not used in the present paper, it is useful to know the form of 0(e 2 ) corrections. The correction terms 
take a complicated form in general, but for the spatially conformally flat metric (7^ = Sij) with (3 l = 0, they take a 
relatively simple form. The components of the Einstein tensor in this case, valid to the accuracy of 0(e 2 ), are given 
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by 



no 



1 2 

"AT 3 ' 



(A25) 



Energy-momentum tensor 



The four- velocity is given by 



u° = [Af 2 -((3 k +v k ){p k +v k )\ 
vl 



feM-i/2 



uV. 



where v l is the spatial velocity, and 



u Q = -u°[N 2 - (3 k {p k + v k )} , 
Ui = u" (vi + Pi) , 



and Vi = ~/ijV J . The components of the energy momentum tensor are then given by 

T° = -(u°) 2 (p + P)[M 2 -(3 k (v k +p k )]+P, 
T°i = (u°f(p + P)(v l +f3 t ) , 

= (u°) 2 (p + P)v i (v j + l3 j ) + 6 i j P. 



(A26) 
(A27) 

(A28) 
(A29) 



(A30) 
(A31) 
(A32) 



We note that T% = —p + 0(e 2 ) if f3 l and v % are both of 0(e). If we choose our spatial coordinates to be comoving 
with the fluid, we have v l = 0. Then the (0, i)-component of the Einstein equations tells us that T ( \ = 0(e), which 
implies (5 l = 0{e). Hence the (0, 0)-component of the Einstein equations gives a local Friedmann equation at each 
spatial point x l . In other words, as long as we are concerned with Einstein gravity, it is unnecessary to assume 
ft 1 = 0(e), but the only condition we need to obtain the local Friedmann equation is the comoving condition v % = 
for the spatial coordinates. 

Defining the projection tensor 

= g**" + u»u v , (A33) 

i.e., projecting orthogonally to the velocity u M , we can write down the momentum conservation equation 

h Xv V»T^ = , (A34) 

which gives in components in the comoving gauge ((3 l = v l = 0) 

(p + P) Di \nN + DiP = , (A35) 

which reduces to the linear result in the comoving gauge on large scales. It may worth noting that this holds for 
general jy. 

The Raychauduri equation [43| in the comoving gauge ((3 1 — v l — 0) and assuming zero vorticity is given by 



^9+ l -e 2 + l -{ P +'iP) + 0{e 2 ) = Q. 



(A36) 



With the identification of 9 = 3H, this is equal to the time derivative of the local Friedmann equation, which can be 
also obtained by combining the (0, 0)-component of the Einstein equations as given by Eqs. (|A22ll and (|A30|I . and the 
trace of the (i, j)-component as given by Eqs. (|A17|) and i|A32(l . 
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APPENDIX B: EVOLUTION EQUATIONS FOR 7y 

Here we write down the traceless part of the spatial components of the Einstein equations, that is, the evolution 
equations for 7^. They can be found, for example, in Eqs. (2.11) and (2.12) of Shibata and Sasaki (SS) ^||. Their 
Aij is related to our Ay by Ay = e~ 2a Aij, or A l j = j^Aij — A l j. Equation (2.11) of SS is 

dt% = -2AfAij + £fyj - p iJ d k k , (Bl) 
where £p is the Lie derivative along f3 k , given for a second-rank tensor Qij by 

£ Qij = P k d k Q l3 + Q ik d J p k + Q k AP k . (B2) 

Thus, with the assumption that f3 k = O(e), the assumption dciij — 0(e) is equivalent to A^ = 0(e). 
Equation (2.12) of SS is 

dtAtj = N (KAu - 2A lk A k 3 ) + ^ [AT ( (3) i? y - ^ (3) i?) - (DiDjff- ^-D k D k N 
2 -A l3 d k p k -^G^ (.S, ; 



+£pA ij - -Aijdkp* - SttG^ ( Stj - ^-S k k ) , (B3) 



where Sij is the spatial projection of the energy-momentum tensor, 

Sa = T V} , S k k = 7 W 5 M . (B4) 
For a perfect fluid, or in the absence of anisotropic stress, we have 

S tj = (p + P)u lU] + jij P . (B5) 

Hence 



Sij - 1 fS k k = {p + P)(u°f (Vi + + Pj) - -^(v k + &)(«* + /3 fe ) , (B6) 



3 



3 



which is of second order in e. 

Assuming the anisotropic stress is negligible, then Eq. (|B3|) reduces to 

dtAij = UK Aij + 0(e 2 ) = -3 d t a A t] + 0(e 2 ) . (B7) 

Therefore, if A.^ = 0(e), Eq. I|B7(1 has a decaying solution, 

Aij = e- 3a Cij , (B8) 

where Cy = 0(e) and dtCij = 0(e 2 ). Assuming that this decaying solution is absent (or ignorable), we have 
dtjij = 0(e 2 ), which we assumed in the text. 
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